Differentiating the Method of Conjugate Gradients
نویسندگان
چکیده
The method of conjugate gradients (CG) is widely used for the iterative solution of large sparse systems of equations Ax = b, where A ∈ R is symmetric positive definite. Let xk denote the k–th iterate of CG. This is a nonlinear differentiable function of b. In this paper we obtain expressions for Jk, the Jacobian matrix of xk with respect to b. We use these expressions to obtain bounds on ∥Jk∥2, the spectral norm condition number of xk, and discuss algorithms to compute or estimate Jkv and J T k v for a given vector v.
منابع مشابه
Conjugate gradient neural network in prediction of clay behavior and parameters sensitivities
The use of artificial neural networks has increased in many areas of engineering. In particular, this method has been applied to many geotechnical engineering problems and demonstrated some degree of success. A review of the literature reveals that it has been used successfully in modeling soil behavior, site characterization, earth retaining structures, settlement of structures, slope stabilit...
متن کاملLanczos and Linear Systems Lanczos and Linear Systems
abstract Lanczos's major contributions to the numerical solution of linear equations are contained in two papers: \An Iteration Method for the Solution of the Eigenvalue Problem of Linear Diierential and Integral Operators" and \Solutions of Linear Equations by Minimized Iterations ," the second of which contains the method of conjugate gradients. In this note we retrace Lanczos's journey from ...
متن کاملNatural Gradients in Practice: Non-Conjugate Variational Inference in Gaussian Process Models
The natural gradient method has been used effectively in conjugate Gaussian process models, but the non-conjugate case has been largely unexplored. We examine how natural gradients can be used in non-conjugate stochastic settings, together with hyperparameter learning. We conclude that the natural gradient can significantly improve performance in terms of wall-clock time. For illconditioned pos...
متن کاملParallel Preconditioned Conjugate-Gradients Methods on Transputer Networks
We show how highly efficient parallel implementations of basic linear algebra routines may be used as building blocks to implement efficient higher level algorithms. We discuss the solution of systems of linear equations using a preconditioned Conjugate-Gradients iterative method on a network of transputers. Results are presented for the solution of both dense and sparse systems; the latter bei...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- SIAM J. Matrix Analysis Applications
دوره 35 شماره
صفحات -
تاریخ انتشار 2014